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We investigate theoretically 2D electronic transport in Si:P 5-doped layers limited by charged- 
dopant scattering. Since the carrier density is approximately equal to the dopant impurity density, 
the density dependent transport shows qualitatively different behavior from that of the well-studied 
2D Si-MOSFETs where the carrier density is independent of the impurity density. We find that 
the density dependent mobility of the Si:P system shows non-monotonic behavior which is exactly 
opposite of the non-monotonicity observed in Si-MOSFETs — in the Si:P system the mobility first 
decreases with increasing density and then it increases slowly with increasing density above a typical 
density 10 14 cm -2 (in contrast to Si MOSFETs where the mobility typically increases with density 
first and then slowly decreases at high density as surface roughness scattering dominates). In the 
low density limit (or strong screening limit) mobility decreases inversely with increasing density, but 
in the high density limit (or weak screening limit) it slowly increases due to the finite width effects 
of the 2D layer. In the intermediate density regime (1/a < 2k p < qrF, where a, kp, and qtf are 
the confinement width of the 5-layer, Fermi wave vector, and Thomas-Fermi screening wave vector, 
respectively) the density dependent mobility is approximately a constant at the minimum value. 
However, the calculated mean free path increase monotonically with density. We also compare the 
transport scattering time relevant to the mobility and the single particle relaxation time relevant 
the quantum level broadening, finding that the transport scattering time could be much larger than 
the single-particle scattering time unlike in Si MOSFETs where they are approximately equal. 



I. INTRODUCTION 

In semiconductor based quantum computers, where 
electron spin (or charge) localized at dopant sites could 
be used as qubits 1 -, it is extremely important to very 
precisely place dopants with atomic precision within the 
host lattice^. This has motivated a materials science 
initiative to create precisely located dopant structures 
within semiconductor hosts, e.g. P in Si or Ge. In 
particular, the scanning tunneling microscopy (STM) 
based lithography 3 - has been used to position dopants 
with atomic scale precision^ within semiconductors. The 
ability of controlling the location of individual dopant 
atoms within a semiconductor provides future poten- 
tial towards atomic scale devices 1 -^. Recently, Si:P 6- 
doped two dimensional (2D) structures in silicon have at- 
tracted interest with their potential applications in multi- 
qubit quantum computer architectures 2 - and atomic-scale 
device^"— . The effect of doping density on electronic 
transport in Si:P (5-doped layers grown by phosphinc dos- 
ing and low temperature molecular beam epitaxy have 
been investigated in experiments^"— . Unlike other 2D 
systems, the carrier concentration in this (5-dopcd layers 
can reach 10 14 cm -2 or higher, which exceeds a typi- 
cal carrier density found in GaAs-based systems or Si- 
MOSFETs by more than two orders of magnitude. 

Even though the carrier density of Si:P (5-doped lay- 
ers is extremely high the measured mobilities^"— at 4K 
are very low n ~ 30 — 200 cm 2 /Vs in P doping levels 
of ~ 10 14 cm" 2 . Thus, in these 6 layers the mean free 
path / is of the order of only a few nm and the cor- 
responding transport relaxation time r is in the fem- 



tosecond range. In addition, as the carrier density in- 
creases mobility decreases at first, but the mean free 
path increases linearly. These behaviors are opposite 
to those found in extensively studied 2D Si MOSFET 
and GaAs/AlGaAs systems. It is, of course, expected 
that the electronic transport mobility of the Si:P 5-dopcd 
systems is lower than that in modulation-doped GaAs 
systems because the charged dopant impurities in Si:P 
layers are located inside the layer, but the mobility for 
the Si:P system is substantially lower than Si MOSFET 
mobility also. However, the density dependence of low- 
temperature transport in Si:P layers is unusual because 
it is dominated entirely by the dopant charged Coulomb 
impurity scatterin g 13 i 14 . The mobility increases with in- 
creasing density when the charged impurity scattering 
dominates in the high density regime. In Si-MOSFETs, 
due to the surface roughness scattering at high densities, 
the mobility decreases with increasing carrier density in 
the high density regime dominated by surface roughness 
scattering which however plays no role in the transport 
properties of Si:P delta-layers since the carriers are lo- 
cated away from the surface inside the semiconductor. It 
is also reported that in Si:P (5-dopcd samples with the 
doping density of < 10 13 cm" 2 a non-Ohmic behavior 
is observed, suggesting that a transition to Anderson lo- 
calization for Si:P (5-doped layers occurs below n ~ 10 13 
cm" 2 . The similar low-density localization behavior is 
also observed in other 5-doped samples like Si:Sb and 
Si:B (5-doped layers 1 ^. 

These experimental features of (5-doped layers provide 
a unique opportunity to study the carrier density depen- 
dent transport of ultrahigh density 2D systems limited 
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by only charged Coulomb disorder. Since the carriers re- 
main confined to a dopant plane by strong-confining elec- 
trical fields, charged Coulomb scattering by the ionized 
dopants themselves significantly limits the mobility. In 
spite of there being very interesting experimental studies 
in the ultrahigh density <5-doped systems, the systematic 
theoretical investigation of transport properties has not 
yet been carried out. The purpose of this paper is to the- 
oretically study the transport properties of Si:P <5-dopcd 
layers. It is obviously important to understand the de- 
tailed transport properties for eventual applications to 
atomic-scale devices and semiconductor multiqubit ar- 
chitectures. 

In addition to atomic scale device application motiva- 
tion there are several fundamental reasons which make 
transport studies of 2D (5-doped structures interesting 
and important in its own right. First, these systems al- 
low the study of 2D carrier transport properties at un- 
precedented high carrier densities of 10 14 cm" 2 or above, 
which is much higher than the typical 10 10 — 10 12 cm~ 2 
density regime studied extensively in the context of GaAs 
heterostructures and Si MOSFETs. Second, these sys- 
tems allow for studying carrier transport limited almost 
entirely by long-range Coulomb disorder (arising from 
the ionized dopants producing the carriers in the sys- 
tem) without complications from other mechanisms often 
operational in 2D semiconductor systems such as sur- 
face roughness scattering and alloy disorder scattering. 
Third, the (5-doped 2D systems enable the study of car- 
rier transport in the scattering regime where the carrier 
density is roughly equal to (or perhaps even somewhat 
less than) the impurity density since each electron must 
leave behind an ionized impurity. Fourth, by changing 
carrier density, one could study the strong screening to 
the weak screening 2D transport regime, which is not 
possible in other 2D systems. All of these features are 
unique to the (5-doped 2D systems making them com- 
plementary to the traditional 2D semiconductor systems 
and thus providing particular fundamental impetus to 
our work. We mention that 2D electron systems possess 
the highly peculiar behavior, arising from the constant 
2D density of states which is energy independent. The 
low-carrier density regime is the strong-screening regime 
(where qxF > &f) and the high-density regime is the 
weak-screening regime (where qtf < ^f), and therefore, 
the ability to go to very high carrier density provides a 
completely new perspective in the study of 2D semicon- 
ductor systems, which is not available in the well-known 
Si MOSFET and GaAs heterostructure based 2D sys- 
tems. 

In this paper we investigate the low-temperature elec- 
tronic transport in Si:P (5-doped ultra-high density 2D 
electron systems. Employing Boltzmann transport the- 
ory, we calculate density dependent mobility and scat- 
tering times of Si:P (5-doped layers in a variety of ex- 
perimental situations. The electrons in realistic (5-dopcd 
systems, where there are likely sources of additional scat- 
tering arising from unintentional background impurities, 



may localize at low densities (n < 10 13 cm -2 ), but in this 
paper we neglect localization effects and show our calcu- 
lated results in the wide range of density 10 11 < n < 10 15 
cm~ 2 . In fact, our theoretical results based on the semi- 
classical Boltzmann theory give a clue about the den- 
sity range of strong localization in the system, which we 
discuss later in the paper. All our explicit results (as 
presented in our figures) for mobility, mean free path, 
and scattering time explicitly assume that the 2D carrier 
density is precisely equal to the 2D quenched charged 
impurity density, and therefore, our mobility results are 
an upper bound to the experimental mobility since the 
real systems studied in the laboratory are likely to have 
additional impurities not related to the dopants used in 
creating the electron gas. Wc expect the measured mo- 
bility and conductivity to be somewhat smaller than that 
calculated by us, but we expect the experimental mobility 
to approach our theoretically calculated results with im- 
provement in the sample quality. The calculation within 
the density functional theory shows that the band struc- 
ture of a Si:P (5-doped layer is a symmetric V-shapcd po- 
tential welliSr^S, but we use the square well model with 
a width a as a (5-doped 2D layer. We confirm that our 
calculated results are not sensitive to the shape of the 
confinement potential by also doing a calculation using a 
triangular confinement potential. The calculated trans- 
port results depend only on the characteristic confine- 
ment length V defining the thickness of the 2D electron 
layer, and not on the details of the confining wavefunc- 
tion itself. 

The rest of this paper is organized as follows: in sec- 
tion II wc provide the detailed transport theory and the 
analytic results for both high and low densities, and in 
section III we show the numerical results of mobility and 
scattering times of Si:P (5-doped layers. We conclude in 
section IV with a discussion. 



II. THEORY 

The Si:P (5-doped layer is a conducting layer of atomic 
thickness, which is formed in a pure semiconducting crys- 
tal matrix due to dopant atoms located within one crys- 
tallographic plane. The electrons in the potential well 
form a 2D electron gas in the plane of the 8 layer, and 
behave as 2D free carriers with appropriate effective mass 
and spin/ valley degeneracy quantum mechanically con- 
fined by the ionized dopant-induced electric field within 
a quasi-2D layer of thickness V. The electrons in a Si:P 
(5-dopcd layer are spatially confined along the normal di- 
rection with respect to the layer of P dopants. In this 
paper we designate the 2D layer as (x — y) plane and the 
normal direction of the layer as the z direction (z = 
being the center of the well). The wave function of an 
electron in the layer is given by 

V>(r, z) = exp(ik • r)<fi(z), (1) 
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where 4>(z) is the confinement wave function in z- 
direction, k is the 2D electron wave vector in the plane, 
and r = (x,y). The calculated band structure of Si:P 5- 
doped layer depends sensitively on the disorder model 
(i.e., the nature of symmetry, the exact placement of 
dopants, other unknown fixed charges in the system, and 
the size of the unit cell, etc.). However, for simplicity, 
we assume that the confinement profile is described by a 
square quantum well structure. Although the quantum 
well may be a simplification of the actual situation^—, 
our calculated results describe quantitatively the trans- 
port properties of the Si:P (5-doped layer since the details 
of the confinement are not quantitatively important ex- 
cept for the confinement width, which we parametrize. 

To calculate the 2D mobility of Si:P (5-doped layers, 
we use the Drudc-Boltzmann scmiclassical theory for 2D 
transport limited by the scattering from charged impu- 
rities (i.e., dopants) 13 ' 14 . We assume that the impurities 
(P dopants) are randomly distributed in the x-y plane at 
z = 0, i.e., the center of the quantum well and the 2D 
carrier conductivity is entirely limited by charged im- 
purity scattering. Since the Fermi temperature is very 
high (T F ~ 700 K) at the density n = 10 14 cm" 2 and 
the Bloch-Griineisen temperature at this density is very 
high we neglect phonon scattering in this calculation. 
The Bloch-Griineisen temperature at the electron den- 
sity n = 10 14 cm~ 2 becomes Tbg = 2kphv p \ l ~ 250K 
with the phonon velocity of Si v p h = 9.13 x 10 5 cm/s. 
Below Tbg acoustic phonons cease to be appreciably ex- 
cited, and no longer contribute to the relaxation rate 2 ^. 
We therefore include scattering by the random quenched 
charged dopants as the only resistive mechanism in our 
theory. 

By considering the charged impurity centers located 
at z = with an impurity density and taking k and 
k' to denote the 2D electron wave vector before and af- 
ter scattering, respectively, by a Coulomb scatterer (i.e., 
P dopant) we have the transport scattering time in the 
Born approximation 
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where Vi(q) is the screened Coulomb potential for 
electron-charged impurity interaction, q = k — k' is the 
momentum transfer, #kk' is the scattering angle, and 
€k = (Kk) 2 /2m is the electron energy. 

The unscreened potential of a charged center located 

, where K is 



I' 



at (r i , z = 0) is given by Ui (r) = 

the background dielectric function of Si. The 2D Fourier 
transform of Ui(r) becomes Ui(q) = v(q)Fi(q), where 
v(q) = 2ne 2 1 ' nq is the 2D bare Coulomb interaction and 
Fi(q) is the form factor of electron- impurity interaction 
and given by 



Fi(q) = J dze~ 
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The form factor Fi (q) becomes unity in the limit of van- 



ishing q. The screening effect can be included by dividing 
Ui(q) by the RPA dielectric screening function, e(q), due 
to the 2D electrons themselves, 



e(q) = l + v(q)F(q)U(q), 
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where v(q) is the 2D bare Coulomb interaction, F(q) 
is the form factor for electron-electron interaction, 
and n(q) is the 2D finite wave vector polarizability 
functio n 13 ' 14 . Thus, the screened Coulomb potential for 
electron-charged impurity interaction becomes Vi(q) = 
Ui(q)/e(q). The form factor F(q) for electron-electron 
interaction associated with the confinement wave func- 
tion is defined by 
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where 4>(z) is a confining wave function in the z direction. 
When the (5-doped impurities are located at z from the 
center of the confinement potential the electron-impurity 
Coulomb interaction has the factor e -9 ' 2 ' and becomes 
Vi(q,z)=Vi(q)e-*\*\. 

If we take the confining potential in the z direction to 
be of square-well form, then the quantized ground state 

wave function is given by <fi(z) = yf cos(7r ~z / 'a) for \z\ < 

a/2 and 4>(z) = for \z\ > a/2 with a denoting the width 
of the square well. In this case the form factors Fi{q) and 
F(q) are calculated as 
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As a result of the external dc electric field (E) driv- 
ing the electric current in the system the Fermi distri- 
bution function is deformed from equilibrium distribu- 
tion f (e k ) = l/[e~^-^) + 1] where /3 = l/k B T is 
the inverse temperature and no is the chemical poten- 
tial, i.e., /(ejfe) = /o(efe) + g{^k) where g(e k ) is propor- 
tional to the applied electric field. In the relaxation time 



approximatio: 
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-r(e fe )eE • v k df {e k )/de k , 



where v k = Kk/m is the carrier velocity. Thus, an in- 
duced current density j due to the electric field in a S- 
layer is given by 
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where g is the total degeneracy including the usual spin 
degeneracy of two plus any possible valley degeneracy 
arising from the bulk band structure of the system. From 
the definition of the drift mobility fi = j / (enE) and the 
conductivity a = j/E, we have fi = e(r)/m and a = 
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ne 2 (r)/m, where m is the carrier effective mass and the 
energy averaged transport relaxation time (r) is given by 
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where D(e) = gm/(2Trh 2 ) is the density of states with 
total degeneracy g — g s g v , e = (hk) 2 /2m is the 2D elec- 
tron energy dispersion, and r(efc) is the energy dependent 
transport relaxation time given in Eq. ^ . 

Since we are interested in the low-temperature density 
dependent transport (T <C T F ) of Si:P <5-doped layers we 
first consider the zero temperature limit. After perform- 
ing the k' integration in Eq. ([2]) and taking k = kF (i.e. 
the Fermi wave vector) we have 

l -I{qol (10) 
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where qo = Qtf /2k F (otf = g/a>B is a 2D Thomas-Fermi 
wave vector with effective Bohr radius = h 2 n/me 2 ), 
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Due to the form factors it is not possible to get the full an- 
alytic formula for Eq. (|12l) . However, we can get asymp- 
totic forms of Eq. (|12[) for different density regions where 
large or small qo approximation may apply. For Si(100) 
we have ~ 30 A for electrons in the conduction 
band and qo w 4.7g 3 / 2 /y/n~, where h = n/(10 10 cm~ 2 ). 
Thus, in the low density limit (n < 10 12 cm~ 2 ) g « 1 
(i.e., strong screening limit), and if 2k F a -C 1 (i.e. 
a <C l/2k F ), then the form factors are close to unity 
and I(qo) is given by 



(13) 



Now the scattering time becomes r(k F ) ~ 2to<Jo/tt oc 
1/rii, which is independent of carrier density n. However 
in Si:P <5-doped layer the dopant density rii is equal to 
the carrier density n because every P atom gives up one 
free electron to the layer (as stated above, in real sam- 
ples it is more likely that m > n since there could be 
background unintentional impurities in the system, but 
n = rii provides a lower bound for the impurity density). 
Thus, the scattering time t oc l/n and the density de- 
pendent mobility is inversely proportional to the density, 
/x oc l/n, in the strong screening limit (qo 3> 1). 

In the high density limit (n > 10 14 cm -2 which is the 
relevant density range in the currently available Si:P 5- 
doped layers) qo S> 1 (i.e. weak screening limit) we can 
approximately express Eq. (|12p as 
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Then, we have 



tt/2 for 2k F a < 1, 
2/k F a for 2k F a > 1. 
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Thus for 2k F a < 1 we have r(k F ) = tqtt/2 oc n/rij, 
and mobility u ~ n/rii. For rii = n the mobility is 
independent of carrier density. For 2k F a > 1 we have 
r(k F ) = Tok F a/2 oc n?/ 2 /rii, and the mobility becomes 
[i r-j r?l 2 1 rii. For m = n the mobility increases as y/n as 
the carrier density increases. 

Combining all the results we have so far we can con- 
clude that as the carrier density (or doping density) in- 
creases the mobility decreases as l/n up to 2k F ~ qr F - 
As the carrier density increases further beyond this point 
the density dependent mobility saturates up to 2k F ~ 
1/a, and then it increases as \fn when 2k F > 1/a. Thus 
the mobility behavior in the delta-doped layers as a func- 
tion of carrier density (or dopant density since we assume 
them to be the same in our model) has three distinct 
regimes: a low-density strong-screening regime where the 
mobility decreases (as l/n) with increasing density, an 
intermediate-density saturation regime where the mobil- 
ity is approximately a constant in density, and then a 
high-density weak-screening regime where the mobility 
rises with density (as ji 1 / 2 ). This behavior, which is 
clearly manifested in our full numerical results shown 
in the next section, is in sharp contrast with 2D trans- 
port in Si MOSFETs where the mobility first increases 
with carrier density as Coulomb disorder gets screened 
out and then decreases with carrier density as surface 
roughness scattering becomes dominant. We note that 
if the system is strictly two-dimensional (i.e. the elec- 
trons are confined in an infinitesimally thin 2D layer), 
then the form-factors in Eq. (fT2"j) both become unity, and 
the mobility will decrease monotonically with increasing 
carrier density finally saturating at high enough density 
where go "C 1 condition is satisfied. The third regime 
of an eventual n 1 / 2 increase of mobility at high density 
arises entirely from the quasi-2D form factor effect which 
becomes crucial at very high density when k F 3> 1/a, 
which would not apply if a = 0, i.e. in the strict 2D limit. 
We emphasize again that our theory is an upper bound 
to the expected experimental mobility since we are only 
including scattering by the charged dopants themselves 
ignoring all other possible scattering sources. 



III. RESULTS 

In this calculation we use the following parameters; 
for Si(100) surface orientation the effective mass m* = 
0.19m e and valley degeneracy g v = 2, and for Si(110) 
surface orientation m* = 0.4m e and g v = 4. We use the 
spin-degeneracy g s = 2 everywhere. We set the carrier 
density to be equal to the dopant impurity density, i.e., 
n = n,. We also take a square quantum well with width 
a as a confinement potential. (We have explicitly verified 
that our numerical results change little for a triangular 
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FIG. 1. Calculated mobility of Si:P 5-doped layers as a 
function of density at T = for different quantum well width 
a — 10, 5, and 2nm. Black lines indicate the results calculated 
with Howard-Fang (HF) variational wave function. Here im- 
purities are located at z — (i.e. the center of quantum well). 




FIG. 2. (a) The calculated mobility in a Si:P <5-doped layer 
as a function of density at a finite temperatures, T — 5K. 
The black (red) lines represent the mobility calculated with 
parameters of Si(100) [Si(110)]. The well width a = 2 nm is 
used in this calculation, (b) The exponent of density depen- 
dent mobility calculated in (a), i.e. exponent in the relation 
of n oc n a . (c) fei as a function of density where the mean 
free path is given by I = vft. 



potential well as long as the same effective quasi-2D con- 
finement width is used in the theory.) We also mention 
that mobility is simply related to the conductivity a by 
the formula a = nefj.. 

In Fig. [T] we show the calculated mobility of Si:P 6- 
doped layers as a function of density for different quan- 
tum well width a = 2, 5, and 10 nm. For comparison 
we also show the results (black lines) calculated with 
Howard- Fang (HF) variational wave function 1 ^. In Fig. Q] 
charged impurities (dopants) are located at z = (i.e. 
the center of quantum well). In the low density limit 
(or strong screening limit, go = qpF^kp 3> 1) the cal- 
culated mobility is inversely proportional to the density 



(fx ~ V") as the density increases upto n\ ~ gq^p/lftir. 
In this limit the effect of the finite 2D layer width is 
small since kpa is small and the form factors are effec- 
tively unity. Note that for Si(100) surface orientation 
q TF w 0.13 A" 1 , which gives n x ~ 1.3 x 10 13 cm" 2 . 
For Si(110) surface orientation qpF ~ 0.52 A -1 and 
m ~ 4.2 x 10 14 cm" 2 . When the density increases fur- 
ther, n > n±, mobility has a minimum value at qo ~ 1. In 
the weak screening high-density limit (go "C 1) the form 
factor of electron-impurity interaction arising from the 
finite width of the 2D layer affects the density dependent 
mobility since a > 1 now and the form factors become 
relevant. For both n > n\ and n S> g/(\6iva 2 ) the mobil- 
ity increases as n 1 / 2 . In Fig.[T]the non-monotonic density 
dependent mobility as a function of density is shown, 
i.e., as the carrier density increases mobility decreases 
first, reached a minimum, and then it increases again. 
Thus, the three density regimes (low, intermediate, and 
high) discussed analytically in the last section are clearly 
apparent in Fig. Q] with mobility decreasing first with 
increasing density, then saturating around a minimum 
value around an intermediate density before increasing 
again at higher density. The characteristic intermediate 
density for the mobility minimum increases as the quan- 
tum well width decreases. Fig. Q] also shows the degener- 
acy dependent mobility. At low densities n < 10 13 cm" 2 
the mobility of Si(110), which has a higher degeneracy 
factor and hence stronger screening, is higher than that 
of Si(100), but at high densities the mobility of Si(100) 
is higher. We can understand this degeneracy depen- 
dence from the scattering time. From Eq. (fTUj) we find^ 2 - 
/i oc g 2 for (jo > 1 (strong screening limit) and /i oc g" 1 
for go -C 1 (weak screening limit). 

In order to clearly see the high density behavior of mo- 
bility we show in Fig. HJa) the calculated mobility in lin- 
ear scale as a function of density at a finite temperature 
T = 5K. The black (red) lines represent the mobility 
calculated with parameters of Si(100) [Si(110)] . The well 
width a = 2 nm is used in this calculation and the impu- 
rities are located at z = 0. Fig.^b) shows the exponent 
of density dependent mobility for the results shown in 
(a), i.e. exponent in the relation of oc n a . As described 
in Sec. II the exponent a of Si(100) samples eventually 
approaches 1/2 as the carrier density increases. Fig.^c) 
shows the calculated kpl as a function of density where 
the mean free path is defined by I = VpT. Note that 
kpl = (2cr / g)(h/e 2 ), where the conductivity a has a unit 
e 2 jh and g is the total degeneracy. In terms of mobility 
we can express it as kpl = 0.83^m/g, where Li is mea- 



sured in units of cm /Vs and n = n/(10 



14, 



2 ). The 



calculated kpl increases super- linearly due to the finite 
width effect. In Sec. II we show that /i oc n 1 / 2 for go <C 1 
and kpa 3> 1. Thus in this high density limit we expect 
that the mean free path increases linearly with density, 
I oc n, and this behavior is experimentally observed in 
Si:P 5-doped layers 1 ^. However, in the strict 2D limit 
kpl increases linearly with density and / oc ^/n. The lin- 
ear increase of I with n arises from the finite width effect 
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FIG. 3. (a) The transport scattering time (black lines), Tt, 
and single particle relaxation time (red lines), t s , as a function 
of density for Si(100) (dashed lines) and Si(110) (solid lines) at 
T = 0. The charged impurities are located at z — (center of 
the well). The well width a = 2 nm is used in this calculation. 
In (b) we show the scattering times assuming the impurities 
be outside the quantum well z = 5nm. 



of the layer as discussed analytically in Sec. II. 

So far we have only considered the transport scattering 
time (t) which is relevant to the conductivity or mobility 
as shown in Sec. II. In general, in the presence of disorder 
scattering there are two distinct relaxation times2^: the 
scattering lifetime or the transport relaxation time (r or 
Tt) and the quantum lifetime or the single particle relax- 
ation time (r s or T q ). These two characteristic times of 
the system differ by the important (1 — cos 9) vertex cor- 
rection factor. The relaxation rate t s _1 is given by mak- 
ing the replacement (1 — cos 9) — > 1 in the integrand for 
the formula for t _1 in Eq. @. In general, t s determines 
the quantum level broadening, 7 = H/2t s , of the mo- 
mentum eigenstates, and physically, it simply represents 
the time between scattering events between electron and 
impurity. Considering the difference between two scat- 
tering times we calculate the single particle relaxation 
time and compare it with the transport scattering time 
in following two figures. We mention that just as the 
transport scattering time r t determines the conductivity 
or the mobility of the system, the single particle relax- 
ation time t s determines level broadening as measured in 
Shubnikov-dc Haas (SdH) measurements — in fact, the 
Dingle temperature of SdH measurements is given by 7 
defined above using r s . In strong screening systems or 
for short-range disorder scattering (as in 3D metals or 
2D Si MOSFETs), r s ~ r t generally, and one does not 
need to discuss two distinct impurity scattering times. 
This is not true in Si:P <5-layers at high density where 
the system is weakly screened, and the two scattering 
times could differ by a very large factor in the presence 
of long-range Coulomb disorder. 

In Fig. EUa) we show the transport scattering time, 
t = T t , and single particle relaxation time, r s , of Si:P 
(5-doped layers as a function of density for both Si(100) 
and Si(110). The charged impurities are located at z = 
(center of the well). The well width a = 2 nm is used in 
this calculation. In Fig.^b) we show the scattering times 
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FIG. 4. (a), (b), and (c) show the calculated scattering times 
at T = with parameters of Si(100). (d), (e), and (f) show 
the scattering times at T — with parameters of Si(110). 
The quantum well width a — 2 nm is used in this calculation. 
(a),(d) The transport scattering time (black lines), Tt, and 
single particle relaxation time (red lines), t s , as a function 
of density for different impurity locations, z = 0, 1, 2, 5nm 
(bottom to top). (b),(e) The ratio r t /r s as a function of 
density for different impurity locations. (c),(f) The ratio Tt /t s 
as a function of impurity location z for different densities. 



assuming the impurities to be outside the quantum well, 
i.e. impurities at z = 5nm. By replacing (1— cos9) with 1 
in Eq. @ we have t s rs 2ro<7o/ 7r x n ~ 1 f° r Io ^ 1 (strong 
screening limit or low density limit). In the opposite limit 
(qo <C 1) we have r s m tqQq oc n -1 / 2 , assuming n, = n. 
Thus, unlike t, as the carrier density increases the calcu- 
lated single particle relaxation time t s decreases continu- 
ously without any upturn which appears in the transport 
scattering time r [see Fig. Ufa)] ■ When the charged im- 
purities are located outside the quantum well (remote 
doping) the very strong enhancement in the transport 
scattering time r is found since large-angle scattering 
(2/cf scattering) by the remote impurities (which is most 
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FIG. 5. (a) The calculated mobility in a Ge:P 5-doped layer 
as a function of density at a finite temperatures, T = 5K. 
The black (red) lines represent the mobility calculated with 
parameters of Ge(100) [Ge(llO)]. The well width a = 2 nm 
is used in this calculation, (b) kpl as a function of density 
where the mean free path is given by I — vft. 



weighted in r) is strongly suppressed by the separation. 
The single particle relaxation time is also enhanced, but 
the enhancement is much smaller than that of r since all 
angles contribute to t s [see Fig. G^b)] . 

In Fig. 2] we show the calculated scattering times and 
the ration of t to r s for both Si(100) and Si(llO). As 
shown in Figs. 0|a) and (d) the transport scattering 
time, r, shows nonmonotonic behavior and strongly de- 
pends on the location of dopant impurities. However, 
the single particle relaxation time, t s , decreases mono- 
tonically as the density increases, and has weaker de- 
pendence on impurity location. In the strong screening 
limit and z = we have r = t s rj 2tq%/pi cx n^ 1 and, 
therefore, t/t s — > 1. At high densities (go S> 1) and z = 
we have t/t s cx y/n, but due to the factor e~ 2kpZ arising 
from the separation between electrons and impurities the 
ratio is much enhanced as shown in Fig. 2) 

Finally, in Fig. 02a) we show the calculated mobility in 
a Gc:P (5-dopcd layer as a function of density at a finite 
temperature, T = 5K. Recently Ge:P <5-doped layers 
have also attracted a lot of interest due to the higher 
mobility in Ge than in Si^ 4 -. In this calculation we use 
the following parameters: effective mass m* = 0.12m e , 
background dielectric constant re = 16, and the valley 
degeneracy g v = 4 for Gc(100) surface and g v = 2 for 
Ge(llO) surface. We use the well width a = 2 nm and 
assume the dopants are located at z = 0. Fig.[5]Jb) shows 
kpl as a function of density where the mean free path is 
given by I = vpT. We find that the mobility in a Ge:P 
(5-doped layer is enhanced by a factor of 10 compared 
with the mobility in a Si:P (5-doped layer. However, the 
overall density dependent transport behavior of the Ge:P 
is very similar to that of Si:P. Due to the large value 
of the effective Bohr radius in Ge (as ~ 100 A) the 
Qo = ItfI^f in Ge:P is much smaller than that of Si:P. 
Thus we find that after reaching a minimum point in the 
density dependent mobility the upturn density is much 
smaller in a Gc:P <5-dopcd layer than that in a Si:P 6- 
doped layer. 



IV. DISCUSSION AND CONCLUSION 

Before concluding, we first want to discuss two issues 
not discussed in the earlier sections with respect to our 
results. These are the important questions of comparison 
between experiment and our theory and the issue of lo- 
calization, which is of course experimentally relevant, but 
is not explicitly included in our semiclassical Boltzmann 
theory. 

First, in comparing our theory qualitatively to experi- 
ment, we note that our theoretical results are at best an 
upper limit on the experimental mobility since we have 
assumed the ionized dopants to be the only source of re- 
sistive scattering whereas in reality there are likely to be 
other sources of disorder (e.g. unintentional and hence 
unknown impurities) in addition to the charged dopants. 
We therefore expect our results to agree with only the 
best available samples where other sources of scattering 
are presumably suppressed. This is indeed the case. The 
current highest measured mobility of Si:P (100) <5-layer 
comes from the Sandia National Lab^ where the samples 
after annealing reached a mobility at 5K of close to 200 
cm 2 /Vs at a carrier density of n ~ 10 14 cm -2 . According 
to our Fig.[TJ the mobility for the Si:P (100) (5-layer varies 
between 200 and 500 cm 2 /Vs which is consistent with the 
experimental data (and is in fact in agreement with it if 
we use the lower value of the 2D confinement width). 
All other existing experimental mobility data in the lit- 
erature fall below our calculated value (but none above) 
because the currently existing samples are not yet opti- 
mized and still presumably contain substantial amount 
of unintentional quenched impurities. This conclusion is 
consistent with the Sandia finding that the mobility in- 
creased substantially (by more than a factor of 2) upon 
annealing which presumably got rid of some of the un- 
intentional impurities in the system. In any case, our 
results presented in this paper should motivate the ex- 
perimentalists in producing density dependent mobility 
data in Si and Ge 5-doped 2D systems so that a careful 
quantitative comparison can be carried out, which is not 
possible at this stage because of the dearth of detailed 
experimental data. 

Second, we can discuss the onset of strong localization 
within our Boltzmann theory by asking where the Ioffc- 
Reggel-Mott criterion^ for strong localization, kpl = 1, 
starts to get satisfied in our system. This will again pro- 
vide us with an intrinsic limit on localization since we 
assume n = m which would underestimate the actual 
localization density in real samples where n < rii is the 
likely situation because of unintentional impurities. We 
note that our calculated kpl [Fig. H£c)] increases mono- 
tonically with density, and thus the localization condi- 
tion of kpl is a low-density condition, which is consistent 
with experimental data where the localization behavior 
is seen only for n < 10 13 cm~ 2 . In fact, for n = 10 14 
cm" 2 , we get kpl ~ 40 in Si:P (100) 2D system, which 
is far from the strongly localized regime. Using our sim- 
ple formula kpl — 0.83/iri derived in Section III, we see 
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from Fig. [TJa) that for n = 10 13 cm -2 (where our calcu- 
lated fi = 200 cm 2 /Vs), kpl ~ 4 which is rather close to 
the strong localization condition. Therefore, it is indeed 
possible that the actual localization in realistic samples 
(where the mobility is likely to be lower than our calcu- 
lated results assuming that the ionized dopants are the 
only scattering source) occurs around n ~ 10 13 cm" 2 
carrier density. 

We also comment on our finding of a mobility minimum 
at a characteristic density n c in the 2D <5-doped system. 
Since n c defines the crossover from the strong screening 
to the weak screening behavior in the system, we can 
crudely obtain n c from the condition qo = 1, i.e. qrF = 
2k F . This happens for Si(100) system at n c = 1.2 x 10 13 
cm" 2 , which is roughly where the mobility starts to go 
from a decreasing function of density to an increasing 
function of density in Fig. [T] It will be very interesting 
to observe this strong to weak screening crossover directly 
exp er iment ally. 

In conclusion, we have studied the effect of doping den- 
sity on electronic transport in 2D Si:P <5-doped layers. 
Since the carrier density is equal to the dopant impurity 
density we find that the density dependent transport is 
qualitatively different from that of Si-MOSFETs where 
the carrier density is independent of the impurity density 
since carriers are induced by a gate and the impurities are 
unintentional oxide charges at the interface. The density 
dependent mobility of a Si:P (5-doped layer is also differ- 
ent from that of a modulation <5-doped 2D GaAs quantum 
well system because of the physical separation between 
electrons and dopants in the GaAs system. We find that 
the density dependent mobility of the Si:P system shows 
a non-monotonic behavior. At qo = qxF/^kp ~ 1 the 
mobility has a minimum value. However, the calculated 
mean free path increases monotonically with increasing 
density as observed experimentally-^. We also calculate 
both the transport scattering time r and the single par- 
ticle relaxation time t s in Si:P <5-doped layers. The r 
shows nonmonotonic behavior like mobility and strongly 
depends on the location of impurities, but the t s keeps 
decreasing with increasing density and is less sensitive to 
the location of impurity. 



In this paper we have simplified the confinement po- 
tential as a square quantum well and neglected two im- 
portant features of high density Si:P samples; the local- 
ization effects and the multisubband effects. As shown 
in Fig. Q] our results are not sensitive to the shape of the 
confinement potential, and given that the exact impurity 
and dopant configurations in real samples are not known 
accurately, wc do not see much point in trying to improve 
our confinement model. Since the localization effects are 
observed at very low densities (n < 10 13 cm" 2 ) our re- 
sults may not be applicable to the actual experiment in 
this density range. Usually in Si:P i5-doped layers only 
a few subbands are filled. Even though the mobility of 
carriers depends on the intcrsubband scattering^ we ex- 
pect that our results quantitatively explain the observed 
experimental data. In particular, as emphasized through- 
out this paper, our work provides the upper limit on the 
possible achievable mobility in delta-doped 2D systems 
since we include the only source of scattering that must 
always be present in the system, namely, the scattering 
by the ionized dopants contributing the electrons to the 
system. The experimental mobility can be lower than 
our calculated results, but it can never be higher. Our 
estimate of the localization density using the Ioffc-Regge- 
Mott criterion indicates that high-quality intrinsic delta- 
doped samples with very little unintentional background 
impurities may very well show effective metallic behav- 
ior at carrier density 10 13 cm -2 or lower, but much of 
the strongly metallic behavior most likely manifests for 
density around 10 14 cm" 2 or higher. Our work should 
motivate detailed experimental transport studies for this 
interesting and novel class of 2D electron systems which 
allow the study of 2D transport at very high carrier den- 
sity along with very high impurity density, a transport 
regime little studied before in 2D electron systems. 
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